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Charged and neutral B decays into the n + n~K, K + K~K and K ( ^K^K^ systems have been studied. 
The 7T + 7T~ and K + K~ final state interactions in the 5-wave isoscalar state are treated in an unitary 
way. The long-distance contributions called charming penguins are introduced in the model. Ef- 
fective mass distributions, branching ratios and some asymmetries are succesfully compared with 
the recent BaBar and Belle data. A particularly large negative direct CP-violating asymmetry in 
the charged B decays into /o(980)A' is evaluated for one set of the charming penguin amplitudes. 
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1. Introduction 

A good description of final state interactions is crucial to obtain a precise determination of the 
Cabibbo-Kobayashi-Maskawa matrix elements. Recent experimental results on the B decays into 
three pseudoscalar mesons are a rich source of information about weak decay amplitudes and two 
body mesonic strong interactions. Direct CP-violation effects can appear in these reactions. Here 
we study the B decays into titzK and KKK channels. An important part of these processes is a 
production of the %% or KK pairs in the isospin zero S-wave. A very prominent maximum of the 
/o(980) resonance has been seen in the %% spectra obtained by the BaBar and Belle collaborations 
(see, for example refs.[l-4]). One can expect that the other scalar resonance /o(600) plays some 
role at lower %% effective masses. Both resonances are incorporated in a natural way in the unitary 
three-channel model of pion-pion, kaon-kaon and four-pion interactions [||]. We use this model to 
calculate the final state interactions in the different B-decay channels. 

The penguin diagrams are essential in studies of B decays into tltiK and KKK. However, 
different penguin amplitudes interfere destructively in the decay B — ► /o(980),K. This leads to 
much too small branching ratios calculated in the factorization approach. Thus we consider the 
long-distance contributions originating from enhanced charm quark loops. These so-called charm- 
ing penguin terms correspond, for example, to weak decays of B to intermediate DgD^ states 
followed by transitions to the /o (980)7^ final states via cc annihilations. The charming penguins 
have been used in fits of experimental branching fractions for two-body charmless B decays [6-7]. 
Below we show that these long-distance contributions are also needed in analyses of B decays into 
three pseudoscalar mesons titiK or KKK. 

2. Model for the B -> nnK and B -> KKK decays 

We have constructed a model for the following charged and neutral B decays: B^ — > [m^sK^, 
B ± -> (KfysK*, B° -► (kk) s K°, B° -» (KK) S K°, B° -► (kk) s K° and 5° -» (KK) S K°. Here by 
(7l7l)s and (KK)$ we mean n + 7i~ or n Q 7i° and K + K or K°K° pairs in isospin zero S-wave. 

Below we outline only the main features of the amplitudes described in more detail in [Q]. The 
decay amplitudes consist of two parts both related to weak transitions b — ► snn and b — > sss, where 
n denotes u or d quarks. The first part corresponds to the factorization approximation with some 
QCD corrections. The second one is the long-distance amplitude Ald with c- or w-quark loops. Its 
form for the B — > (ti + 71~)sK~ decay reads: 

A LD = [C(m^)rr(m^) +C(m K )r 1 s *{m }tJ[ )} , (2.1) 

where Gp is the Fermi coupling constant, % is a constant which value is close to 30 GeV -1 , mx 
is the kaon mass, T"{m nn ) and T((m KK ) are the non-strange and strange pion scalar form factors 
depending on the effective pion-pion mass m %% . The charming penguin contribution is written 
similarly as in ref . [^] : 

C(m) = - (Ml-m 2 )f„F % (V ub V^P u + V th V*P t ) , (2.2) 

where fnF K = 0.042 GeV, P, and P u are complex pai - ameters multiplied by the products of the 
Cabibbo-Kobayashi-Maskawa matrix elements V. We have performed numerical calculations using 
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two sets of parameters fitted to data on two-body charmless B-decays. We consider the model I in 
which [|] P, = (0.068 ± 0.007) exp[i( 1.32 ±0.10)] and P u = (0.32 ± 0. 14) exp[i( 1.0 ±0.27)], and 
the model II [§] with P, = (0.08 ± 0.02) exp [-/(0.6 ± 0.5)] and P u = 0. 

In our treatment of final state interactions we use not only the pion scalar form factors T"(m nn ) 
and T((m nn ) but also the kaon scalar form factors r 2 "(/n^) and T^m^). The first and the second 
form factors are responsible for the transitions from the nn quark pair and from the ss quark pair to 
the nn pair, respectively. Similarly the third and the fourth form factors correspond to the nn — > KK 
and to the ss — > KK transitions. In the B-decay processes transitions from the nn channel to the 
KK channel and vice versa are possible. They are incorporated in the two-body scattering matrix 
T. The interchannel couplings are included in the following formulae: 

r.">» = /?*•» + £ {k i \^\m)G j {m)T ij {m)\k j ), (2.3) 

7=1 

where \kj) and \kj) represent the wave functions of the two mesons in the momentum space and 
the indices i,j = 1,2 refer to the nn and KK channels, respectively. The center of mass channel 
momenta are given by k\ = y / m 2 /4 — m\ and &2 = yj m 2 /4 — m\. 

In the numerical calculations we use the solution A of the coupled channel model The 
functions Gj(m) are the free particle Green's functions defined in []|] and /?" ' s (m) are the produc- 
tion functions responsible for the initial formation of the meson pairs prior to rescattering. The 
production functions have been derived by MeiBner and Oiler in the one-loop approximation of the 
chiral perturbation theory [§]. 

Let us remark that in this approach a sum of many Breit-Wigner terms typical for the so-called 
isobar model, usually used in phenomenological fits to Dalitz plots, is replaced by a set of the 
coupled meson-meson amplitudes T. These amplitudes can be expressed in terms of the phase 
shifts 8 n7I , 8kk and the inelasticity v\ known from other experiments. Thus no arbitrary phases nor 
relative intensity free parameters for different resonances are needed. In addition, both the T- and 
the T-matrices satisfy unitarity constraints. A generalization of Watson's theorem to multi-channel 
B-decays is thus realized. 

3. Results 

Table 1 shows that the results of our calculations compare reasonably well with the results of 
the experimental analyses averaged by HFAG in spring 2005. 

In Fig. la we compare the model I with the BaBar data [|l|] for the B ± — > ^n'K^ 1 decays. 
One can see that the theoretical solid line describes well the shape of the /o(980) resonance. Fig. 
lb shows the Belle data [Q] in comparison with our model for the neutral B Q decays to n + n~K° 
in a wide range of effective masses. The model provides us with an absolute prediction for the 
branching ratio of the B° decay into 70(980)^° if the effective mass distributions in the range of 
the /o (980) resonance for the charged B decays are reproduced. 

For the model I one obtains a particularly large negative direct CP- violation asymmetry. 
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Table 1: Branching ratios 3§ in units of 10~ 6 , direct CP asymmetry j^cp, time dependent asymmetry param- 
eters s# and 5? for the B— >/o(980)^f,/o(980) — > decay. The model errors come from uncertainties of 
the charming penguin amplitudes. 



decay mode 


observable 


experiment 
HFAG 


Model I 

£ = 33.5 GeV- 1 


Model II 
£ = 23.5 GeV- 1 


B ± 


m 

S^CP 


o 4Q+1.35 
1.26 

-o n+ 019 


8.49 (fit) 
-0.52 ±0.12 


8.46 (fit) 
0.20 ±0.20 


B° 




6.0±1.6 
-0.14 ±0.22 
-0.39 ±0.26 


5.9±1.6 
0.01 ±0.10 
-0.63 ±0.09 


5.8±2.8 
0.0004 ±0.0010 
-0.77 ±0.0004 




0.9 0.95 1 1.05 1.1 



m (GeV) 

7171 v ' 



> 

o 

o 



> 




■10 



- b) 

Belle 


Hi 


... i .... J. ... i . .1 . 


— n 



0.4 0.6 0.8 
m (GeV) 

7171 v ' 



1.2 



Figure 1: n + n mass distributions in the B ± — > n + n K* 1 (a) and in the B° 
lines delimit a band of the /o(980) events. 



k + k K° decays (b). Vertical 
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